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Abstract. A Bayesian analysis of the world's p(y, K + )A data is presented. We adopt a 
Regge-plus-resonance framework featuring consistent couplings for nucleon resonances 
up to spin J = 5/2, and evaluate 2048 model variants considering all possible combina- 
tions of 11 candidate resonances. The best model, labeled RPR-2011, is discussed with 
special emphasis on nucleon resonances in the 1900-MeV mass region. 



1 Introduction 



Dedicated experiments have lead to an impressive increase in the quantity and quality of electromag- 
netic-meson-production data. This situation is likely to improve even further in years to come. With 
such an abundance of experimental information, one anticipates great progress in our understanding 
of the nucleon and its excited states. In order to establish the resonance content of meson-production 
reactions, a number of challenges need to be dealt with: couplings between different reaction channels, 
the occurrence of broad and overlapping resonances, and under-determined data sets that do not allow 
to establish a unique best-fit model. This state of affairs prompts for a further refinement of model 
ingredients, but also calls for a re-evaluation of the machinery used for model selection. 

Electromagnetic production of hyperons (Y) and kaons (K) plays a central role in the search for 
hitherto unidentified resonances that have been predicted by quark models. From an experimentalist's 
perspective, kaon production is appealing, despite its small cross section, because the self-analyzing 
power of hyperons facilitates the measurement of polarization observables. The CLAS collaboration, 
for instance, have recently collected data that constitutes an over-determined experiment (T). 

Because of the relatively high threshold for KY production, the reaction mechanism has contribu- 
tions from a large number of broad and overlapping nucleon resonances (N*). In addition, a crucial 
role is reserved for non-resonant diagrams. This observation constitutes a challenge for single- and 
coupled-channel partial-wave analyses, as well as traditional isobar models; the resonance content of 
N(y, K)Y can only be determined after resonant and non-resonant contributions are cleanly separated. 
The Regge-plus-resonance (RPR) formalism elegantly addresses this issue via a hybrid approach com- 
bining a hadrodynamic description in the resonance region with Regge phenomenology constrained at 
higher energies. 

In this contribution, we discuss a Bayesian analysis of the world's p(y, K + )A data adopting the RPR 
framework. Section|2]summarizes the RPR formalism, and introduces Bayesian-evidence computation 
as a powerful tool for model selection. In Section [5] the optimal resonance content for p(y, K + )A is 
discussed with special emphasis on the 1900-MeV mass region. Finally, we present our conclusions. 
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2 The Regge-plus-resonance framework and Bayesian model selection 

In the RPR approach [2 J, the non-resonant contributions to p(y,K + )A are parametrized in terms of 
K + (494) and K* + (892) Regge-trajectory f-channel exchange. To ensure gauge invariance, the electric 
part of a Reggeized s-channel Born diagram is added. This results in a three-parameter model that 
successfully reproduces data at high energies, and describes the general features of p{y, K + )A in the 
resonance region. 

The model is refined by considering contributions from s-channel N* exchanges using standard 
Feynman diagrams. A key aspect of the current analysis is the use of consistent couplings for high- 
spin (/ > 3/2) resonance exchange. This formalism dictates that each spin- 1/2 resonance-exchange di- 
agram introduces one parameter; high-spin resonances require two parameters. Moreover, on account 
of the adopted multidipole-Gauss form factor 0, a single cut-off parameter suffices in the hadronic 
form factors of the A^-exchange diagrams. As a result, the RPR framework allows for an unambigu- 
ous separation of resonant and non-resonant diagrams while the number of free parameters is kept to 
a bare minimum. The latter is important in the light of Bayesian model selection. 

As mentioned in the previous section, the tools used for model selection are crucial to establish 
the resonance content of meson-production reactions. For a model with a set of parameters a M one 
conventionally defines 



where N is the total number of data points, cr, is the error on data point du and /K^m) is the corre- 
sponding model prediction. The quantity in Eq. [T] obeys a chi-square distribution, and parametrizes 
the likelihood function X-ia/w)- A common approach to model selection involves x 2 minimization 
(equivalent to likelihood maximization) for different model variants, and subsequent comparison of 
the minimal x 2 values. A major drawback of this method is the disregard of Occam's razor. Many 
recipes exist to quantitatively penalize the introduction of superfluous model parameters. They are, 
however, not statistically rigorous. 

An alternative procedure, which we have adopted, is offered by the Bayesian evidence 



where the prior distribution n(aM) expresses any prior knowledge of the parameters' probability dis- 
tribution. The numerical computation of Z is highly nontrivial, and requires a dedicated set of algo- 
rithms |4|. Evaluating the Bayesian evidence for different models M A and M B allows to quantitatively 
express the relative probability of both models given the available experimental data \dk] 15141 



The natural logarithm of the evidence ratio can be interpreted qualitatively with the aid of Jeffreys' 
scale, given in Table I of Ref. |4|. To appropriately handle the statistical and systematic errors from 
various experiments, arescaled Bayesian evidence Z' is introduced 115 141 . 



3 Resonance content of p(y,K + )A in the RPR-2011 model 

In Refs. B5I4I . we presented the details of a Bayesian analysis of the world's p(y, K + )A data in the 
RPR approach as outlined in the previous section. Fig.[T]plots the computed evidence values for 2048 
model variants constructed by combining the 11 candidate nucleon resonances with spin J < 5/2 
listed in Table V of Ref. 0. The model with the highest evidence, referred to as RPR-201 1, has 14 N* 
parameters and features the S n(1535), Sn(1650), Fi 5 (1680), Pi 3 (1720), Pi 3 (1900), F 15 (2000), and 
the missing Di3(1900) and Pn(1900) resonances 0. 
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Za = P(M A \ {d k }) 
Zb P(M b \ {d k }) " 
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Fig. 1. The evidences (- In .3') of the 2048+1 model variants in the RPR model space (blue circles, see text), as a 
function of the number of N* parameters. The best model for a fixed number of parameters is indicated with a red 
square. The RPR-2011 model, which has the highest evidence, is marked with an open triangle. The open circle 
denotes the RPR-201 1 model supplemented with the P lt (l924) resonance. 



While the existence of some resonances, e.g. Sn(1650) and Pi3(1720), is confirmed in various 
reaction channels by many models, the status of resonances in the 1900-MeV mass range is unsettled. 
An overview of considered A^*s in this mass range is given in Table [T] Using Bayesian inference, we 
find decisive evidence for the Pi3(1900) state that was promoted from two- to three-star status in the 
latest PDG update J6). In addition, RPR-2011 boasts two states labeled as "missing". Indeed, we find 
for the Pn(1900) and Di3(1900) decisive and significant to strong evidence respectively ||4). 

The presence of a resonance with J" — 1 /2 + and a mass close to 1900 MeV, might point towards 
the contribution of a bound state of the KKN system first found in Ref. J7). This Pi i (1924) resonance 
was confirmed in other model calculations [8|, and predicted to have a significant decay in the KA- 
production channel [9]. The Bayesian methodology that we adopt allows to quantitatively express the 
conditional probability P(R\D) (P(^R\D)) that a resonance R is (not) needed, given experimental data 
D. Evaluating the probability ratio P(R\D)/P(^R\D) for the Pi i( 1924) resonance given the present 
data would require us to calculate the Bayesian evidence X' of an additional 2048 model variants - a 
CPU-time consuming endeavor. In order to get an impression of a possible role of this resonance, we 
have limited ourselves to the evaluation of X! for the RPR-201 1 set of resonances supplemented with 
the Pn(1924). This model variant encompasses 15 N* parameters, and is marked in Fig. [T] with an 
open circle. Adding the Pn(1924) state to the RPR-2011 set results in the best 15-parameter model, 
but does not surpass the evidence of RPR-2011. We find AlnZ,' = 6.17, and thus decisive evidence 
in favor of the RPR-201 1 model according to Jeffreys' scale. This finding does not exclude a role for 
the Pn(1924) resonance in p(y, K + )A, since an alternative set of resonances including the Pn(1924) 
might outperform RPR-2011. 
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Table 1. The considered nucleon resonances in the 1900 MeV mass range given in the notation Lai2j{M), along 
with the PDG status, spin-parity {J"), Breit-Wigner mass (M) and width (T). 



Resonance 


PDG status 


J" 


M (MeV) 


r(MeV) 


D 13 (1900) 
Pn(1900) 
Pi 3 (1900) 
P„(1924) 


missing 
missing 

missing 


-ill- 
l/l* 
3/2+ 
l/2 + 


1895 
1895 
1900 
1924 


200 
200 
500+ 80 
60 



4 Conclusions 

The RPR framework merges Reggeized f-channel exchange with tree-level s-channel ^'-exchange 
contributions from an isobar approach into an economical model for kaon photoproduction in and 
above the resonance region. The RPR formalism clearly separates non-resonant from resonant ampli- 
tudes which is an asset when searching for the properties of resonances which contribute to p(y, K + )A. 
We have used Bayesian inference to perform model selection. 

Starting from a set of eleven candidate resonances, the optimum model for p(y, K + )A, dubbed 
RPR-2011, is determined. This model encompasses the resonances 5n(1535), Sn(1650), Fi5(1680), 
P B (1720), Pn(1900), Fi 5 (2000), D 13 (1900), and P B (1900). Decisive and significant to strong evi- 
dence is found for two states in the 1900-MeV mass range that are not listed in the PDG. Finally, 
we have explored the possible role for a KKN bound state with J" = l/2 + . Our initial investigation, 
whereby the Pn(1924) resonance is added to the RPR-201 1 set, has resulted in a 15-parameter model 
with the third highest evidence of all 2048+1 considered models. Nonetheless, given the current data 
our analysis indicates decisive evidence in favor of the N* combination of the RPR-201 1 model. 
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